Exploring strange nucleon form factors on the lattice 
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We discuss techniques for evaluating sea quark contributions to hadronic form factors on the 
lattice and apply these to an exploratory calculation of the strange electromagnetic, axial, and scalar 
form factors of the nucleon. We employ the Wilson gauge and fermion actions on an anisotropic 
24 3 x 64 lattice, probing a range of momentum transfer with Q 2 < 1 GeV 2 . The strange electric 
and magnetic form factors, G S E (Q 2 ) and G B M (Q 2 ), are found to be small and consistent with zero 
within the statistics of our calculation. The lattice data favor a small negative value for the strange 
axial form factor G\(Q 2 ) and exhibit a strong signal for the bare strange scalar matrix element 
{iV|ss|iV)o. We discuss the unique systematic uncertainties affecting the latter quantity relative to 
the continuum, as well as prospects for improving future determinations with Wilson-like fermions. 

PACS numbers: 11.15.Ha 12.38.-t 12.38.Gc 14.20.-c 
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I. INTRODUCTION 



The strange quark is the lightest non- valence quark in the nucleon, and as such provides a unique window into the 
structure of the proton and neutron. Lattice QCD represents at present the only first-principles predictive method 
to determine such contributions directly from the underlying theory of the strong interaction. The computational 
framework for doing this is well established with no fundamental barriers to success. Until recently, however, the 
calculation of the required quark-line disconnected diagrams was simply too computationally demanding to provide 
statistically significant results with all uncertainties under control. This is beginning to change. With recent algo- 
rithmic advances and continued increases in available computer resources, a new era is dawning where these effects 
might be determined with precision well beyond both experiment and phcnomcnological estimates. Here we report 
on some recent progress toward this goal. 

Strange nucleon form factors represent an attractive test case because they have also been the subject of a vigorous 
experimental program. In particular, a number of collaborations have sought to measure the strange electric and 
magnetic form factors via parity- violating electron scattering [JQ, notably SAMPLE, A4, HAPPEX, and GO. Recent 
combined analyses 0, H| find values for G S E (Q 2 ) and G S M (Q 2 ) that are small and consistent with zero in the range 
of momenta so far explored. Also of interest is the strange axial form factor G S A (Q 2 ), to which electron scattering 
experiments are relatively insensitive. At present, the best constraints come from the two-decades old neutrino 
scattering data of the E734 experiment at Brookhaven Q . A recent analysis Q , combining these results with those of 
HAPPEX and GO, favors a negative value for G S A (Q 2 ) in the range 0.45 < Q 2 < 1.0 GeV 2 . These may be compared 
with the recent MiniBooNE result Q, which is compatible but at the same time consistent with zero. 

A special case is presented by the strange axial form factor at zero momentum transfer, G S A (0) = As, which may be 
identified with the strange quark contribution to the spin of the nucleon. This quantity is of particular importance, 
given the role sea quarks are thought to play in resolving the proton "spin crisis" In principle, it is accessible in 
deep inelastic scattering, where it is given by the first moment of the helicity-dependent structure function As(x). In 
practice, however, determining the first moment requires an extrapolation of the experimental data to small values 
of Bjorken x, where uncertainties are less under control. There is some tension between the two most recent analyses 
from HERMES 0, which rely on different techniques; the former favors a negative value for As while the latter 
finds a result consistent with zero, within somewhat larger uncertainties. 

Unlike the strange electromagnetic and axial form factors, the strange scalar form factor G S S (Q 2 ) is not directly 
accessible to experiment. At zero momentum transfer, this quantity corresponds to the strange scalar matrix element 
(iVjssjiV}. Often considered in relation to the pion- nucleon sigma term [11] . it is an important parameter in models 
of the nucleon. We also note the pivotal role it plays in the interpretation of dark matter experiments. Many models 
of TeV-scale physics, including the MSSM, yield a dark matter candidate (e.g., neutralino) that scatters from nuclei 
via Higgs exchange. The Higgs is believed to predominantly couple to strange quarks in the nucleon, since the lighter 
quarks have proportionally smaller Yukawa couplings while the heavier quarks are too rare to contribute significantly. 
It follows that the cross section is particularly sensitive to the strange scalar matrix element, which enters through 
the parameter 

m s (N\ss\N) 

/Ts = — m~ n — • (1) 

As emphasized recently in [l^ - [T5l | , Jts is poorly known at present and represents the leading theoretical uncertainty 
in the interpretation of direct detection experiments. A commonly used estimate is that of Nelson and Kaplan (la . Il7j ] 
(via fl4IO. who find fr s = 0.36(14). More recent analyses suggest that the quoted uncertainty may be underesti- 
mated [ljg], and as we discuss further in Section TlV Al recent lattice determinations favor a much smaller value. If 
these results prove to be robust, they would imply that the strange quark is in fact not the dominant contribution, 
and predicted cross-sections should be substantially smaller as a result (l8| . 

In recent years, there has been a great deal of progress in probing the structure of the nucleon on the lattice (l^.l20j. 
With few exceptions, however, such studies have been restricted to the determination of isovector quantities or 
otherwise neglect the contribution of "disconnected diagrams," due to the large cost associated with their evaluation. 
Nucleon matrix elements involving the strange quark are inherently disconnected, making them an excellent test case 
for tackling this challenge. Such a matrix element is shown schematically in Figure [TJ by "disconnected," we mean 
that the diagram includes an insertion on a quark loop that is coupled to the baryon correlator only via the gauge 
field. This requires the calculation of a trace of the quark propagator over spin, color, and spatial indices. Since 
an exact calculation would require a number of inversions proportional to the spatial lattice volume, the trace is 
generally estimated stochastically, which introduces a new source of statistical error whose reduction is discussed in 
Section IITAI 

There has been a resurgence of interest lately in computing disconnected form factors on the lattice, building on the 
pioneering studies of a decade ago [21^ 29] , which were mainly carried out in the quenched approximation. Recent work 
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FIG. 1. Schematic representation of a disconnected diagram, giving a strange form factor of the nucleon. Here T is the 
appropriate gamma insertion for the form factor of interest, and N is an interpolating operator for the nucleon. 



has included a determination of the nucleon's strange electromagnetic form factors [30j . preliminary determinations 
of As [3ll - [3"5l | . and several studies of the strange scalar matrix element [3l|, [33-37]. The latter may also be determined 
indirectly from the quark mass dependence of the nucleon mass via the Feynman-Hellmann theorem (38l-l40| or SU(3) 



chiral perturbation theory [41[. Likewise, a complementary approach for determining the strange electromagnetic 
form factors relies on combining lattice data for connected form factors with chiral perturbation theory using finite 
range regularization [42-44]. In this work, we present a direct determination of the strange form factors on relatively 
large, two- flavor, anisotropic lattices. Some preliminary results were presented in [3lll45j. 

The paper is organized as follows. In Section Hi A[ we discuss general considerations for computing the disconnected 
trace and describe the particular approach employed in our calculation. We discuss our approach for extracting form 
factors from the corresponding matrix elements in Section l!! Bl and give further details of the calculation in Section InTl 
In Section HVl we present our results for G S S (Q 2 ), G S A (Q 2 ), G S E (Q 2 ), and G S M (Q 2 ). We conclude in Section IVl with 
remarks on the steps we are taking to improve future determinations of these quantities. 



II. METHOD 



A. Evaluating the trace 

As discussed in Section U the evaluation of a disconnected form factor requires the trace of the quark propagator, 
times some combination of Dirac gamma matrices (r) over a time-slice of the lattice. The standard method for 
estimating such a trace relies on calculating the inverse of the Dirac operator D against an independent set of random 
vectors r/ Q ),a = 1, • • • , N for each spin (i), color (a), and spatial (x) degree of freedom, 

1 N 

TviTD- 1 ) « - £ n'^TD- 1 ^ , (<nl atMv ) n = 6 xy 6 ab 6 l3 . (2) 

Given a finite ensemble of such noise vectors, this procedure introduces a new source of statistical error, ±a/y/W, 
measured by the variance a 2 . We note that if we choose an appropriate basis, T = (1, 75, 17^, 757^, icr^), for the 
Euclidean gamma matrices, 75 invariance implies all the traces in Eq. [5] are real. Introducing Gaussian random 
vectors, the variance for the real part of 0(rj) for large N is 

° 2 = \ ({[0*(V) + 0( V )} 2 ) n (O'fa) + 0( V )) 2 ) 

= i^tr[i?t-i(. Tj2/ ) jD -i(y ;X )] + I^ tr [r^- 1 ( a; ,y)r^- 1 (y, a :)], (3) 

x,y x,y 

where the operators are 0(rf) — rfTD~ lr q and tr[- • •] now stands for the trace in spin and color alone. It is also 
common in practice to introduce "unitary" noise elements in U(l) or Z2 46], such that rj*rj = 1 which eliminates the 
diagonal terms corresponding to x = y. In either case, the off-diagonal terms for the first term on the right hand side 
of Eq. [3] fall off exponentially in the space-time separation (|a; — y\) dictated by the lowest Goldstone pseudoscalar 
mass in the quark- antiquark channel, giving a divergent variance in the chiral limit. Taking the real part reduced this 
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divergent term by 1/2 and substituted more rapidly falling correlators in the second term, except for the pseudoscalar 
case where r = 75. In this article we employ two modifications to reduce the variance. Additional methods will be 
explored in a subsequent publication. 

First, we choose a random SU(3) gauge transformation Q, x as our unitary stochastic vector, with elements drawn 
from a uniform distribution according to the Haar measure. We also treat the spin contractions exactly, along the 
lines of the "spin explicit method" described in [13]. As a result, Eq. ([2]) is replaced by 

1 N 

Tr(rL)- 1 ) « — tr^^rZ?- 1 ^")] , (0+ a6 0^) n = S xy 5 ad 6 bc /3 , (4) 

a=l 

and Eq. © by 

° 2 = ^E trc ( trs ^ t_1 ^'y) rt ] tT*lrD- 1 (y,x)})+~^tr c (tr s [rD- l (x,y)}tr s [rD- 1 (y,x)]) , (5) 

x=£y x^y 

where tr c and tr s are color and spin traces respectively. Note that all 12 spin/color components in the local term 
with x = y are removed from the variance. In addition, the explicit spin sum results in separating the SU(3) gauge 
trace from the Dirac (spin) trace for each propagator. The variance depends on the individual gamma structure and 
as a result in general falls off faster for large \x — y\ than the lowest Goldstone mode. To determine the specific 
quark/anti-quark channel that contributes, one must perform a Fierz transformation on each term to put the gamma 
matrices in the conventional position for a meson two-point function. There will be only one linear combination 
affected by the Goldstone mode. The other 15 combinations are determined by massive meson channels even in the 
chiral limit. We leave a more detailed analysis to a future publication dealing with the light quark sector where this 
becomes a more critical issue. 

Second, we introduce dilution to reduce the variance, by dividing the stochastic source into subsets and estimating 
the trace on each subset separately [H, [4{|. As a simple illustration consider even/odd dilution, which involves two 
subsets, 

N N 

TriTD- 1 ) « - ]T vi^TD- 1 ^ + - £ ^TD'^ , (6) 

a—l a—1 

where r)e and rji a ' are non-zero only on the even and odd sites, respectively, and rji + — rf 01 ^ gives the original 
noise vector. This may be generalized to a more aggressive dilution pattern where a larger number of diluted sources 
is used, each more sparse. The lighter the quark mass the more aggressively one should use dilution. We combine 
such a scheme with SU(3) unitary noise and an exact treatment of the spin sum. Note that had we instead used 
dilution over the color index, the resulting variance would no longer be gauge invariant. 

A full calculation involves two sources of statistical error: the usual gauge noise and the error in the trace. In this 
investigation for the strange quark, we largely eliminate the latter error by calculating a "nearly exact" trace on each 
of four time-slices with very aggressive dilution. This is accomplished by employing a large number of sources (864 x 12 
for color/spin on a 24 3 x 64 lattice) where each source is nonzero on only 16 sites on each of the four time-slices. The 
sites are chosen such that the smallest spatial separation between them is 6^/3a s . With this aggressive dilution pattern, 
we find that it is sufficient to use a single SU(3) source per subset, provided automatically by the random gauge noise 
in stochastically independent configurations of our ensemble. Any residual contamination, which we observe to be 
small, is gauge-variant and averages to zero. We note that apart from our use of dilution, this approach corresponds 
to the "wall source without gauge fixing" method employed in some of the earliest investigations of disconnected form 
factors [Hill. 

B. Lattice correlation functions 

In Minkowski space, the familiar Dirac and Pauli form factors of the nucleon, Fi(Q 2 ) and -F^Q 2 ), are implicitly 
defined by 



(N(p')\J,\N(p))=u(p') 



7^i(Q 2 ) + ^^ 2 (Q 2 ) 



«(P) • (7) 



Here \N(p)) is a nucleon state with momentum p, u{p) is a nucleon spinor, and we define Q 2 = —q 2 , where q — p' — p 
is the 4-momentum transfer. It is often convenient to consider instead the Sachs electric and magnetic form factors, 
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given by 



G E {Q 2 ) = Fl {Q 2 )-^F 2 {Q 2 ) (8) 



and 

G M (Q 2 )=F 1 (Q 2 )+F 2 (Q 2 ), (9) 

respectively. The contribution of an individual quark flavor (e.g., the strange quark) is defined by replacing the full 
electromagnetic current that appears in Eq. (J7J), 

2 1 - 1 

= gW7 M u - -dy^d - -S7 M s + . . . , (10) 

by J* = s7 M s. The corresponding Sachs electric and magnetic form factors are denoted by G S E (Q 2 ) and G S M (Q 2 ), 
respectively. 

Similarly, the strange quark contribution to the axial form factor of the nucleon, G S A (Q 2 ), is implicitly given by 

(N(p')\s W s\N{p)) = u(p') [ llil5 G s A {Q 2 ) + ^- l5 G P {Q 2 )] u{ P ) . (11) 

In this equation G S P (Q 2 ) denotes the strange quark contribution to the induced pseudoscalar form factor of the 
nucleon, which we will not consider further here. Finally, we note that the strange scalar form factor is trivially given 
by G S S (Q 2 ) — (N(p')\ss\N(p)) . Our main focus will be on the special case of Q 2 = 0, with the matrix element denoted 
simply by (iV|ss|iV). 

Our task is to extract these four quantities from appropriately defined Euclidean correlation functions on the lattice. 
We begin by defining the usual two-point function for the nucleon, with momentum q, 

G (2 \t,t ;q) = (l + j 4 r^e^ s (N^x,t)N a (6,t )). (12) 

X 

Here N a = e a bc(ua Cy^d^u" is the standard interpolating operator for the proton, with smeared quark fields, and 
(1 + 74) projects out the positive-parity state. In everything that follows, we always double our statistics by making 
use of the invariance of the action under time reversal. More concretely, in this case we combine the "forward- 
propagating" correlator G^- 2 \t,to;q) with the backward-propagating G_ (to,t;q), where the "— " subscript indicates 
that (I+74) in Eq. (fT2l has been replaced by (1 — 74). 

Next, we define various three-point functions G^S'(t, t',to;q), where X = S, A,E, M correspond to the disconnected 
scalar, axial, electric, and magnetic form factors, respectively. These are given by 

Gf(t,t',t ;q) = (l + 7 i ) a ^e**{Nfi(iBMW(2'J) ~ (W>@ ,t'))]N a (0,t )) (13) 

x,x' 

for the scalar, 

1 3 

Gf{t,t\t -q) = i^^e i ^'[-i(l+j A ha 5 \ a0 (N^S,t)[A i (x',t')- (A 4 (f,O)]^(0,i )) (14) 

i—1 x,x f 

for the axial, 

Gf(t,t',t Q ;q) = (l+ 74 )^^e i «"- 2 '(iV^(f,t)[V 4 (f ,*')- (V 4 (x> ,t'))]N a (0,t o )) (15) 

x,x' 

for the electric, and 

G$(t,t',t ;q) = -L^e ijfc i^e^'H(l + 74)7*75]°" (N^x,t)[V k (x" ,t') - {V k (x" ,t'))]N a (6,t )) (16) 

^ i,j,k ^ x,x' 

for the magnetic. Here and denote the vector and axial currents, and n q in Eq. (|16p simply counts the number 
of nonzero components of q. For , we utilize the conserved current for the Wilson action, given by 

V„{x + a^fi/2) = - [${x + o M A)(7*. + ^Wtix^ix) + i>(x){^ - l)U^(x)i/;(x + a„tt)] . (17) 
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The lattice spacing carries a label u because we will consider anisotropic lattices for which the temporal lattice 
spacing, 04 = a tl differs from the spatial lattice spacing, a x = a 2 = a 3 = a s . For convenience, we define V^(x) 
on a given site of the lattice by averaging those terms involving the adjacent forward and backward links, V^(x) = 
[Vfj,(x + a^/2/2) + V fJ- (x — a /J /t/2)]/2; since the spatial index determines the phase in the Fourier transform, this 
corresponds to an 0(a 2 q 2 ) redefinition of the three-point function. For the axial form factor, we will present results 
computed both from the analogous point-split current, 

A<-p- s -\x + a^fi/2) = % - [$(x + a^^Ulix^x) + ^(x)^ 5 U„(x)i;(x + a^)] , (18) 

and the standard local current, A^ ocai \x) — i^ix^f/y^ix). 

Note that in Eq. (Tl^| - p^|) we always employ the vacuum-subtracted value of the current, [J(x,t) — (J(x,t))], even 
though this is only strictly necessary when J is the scalar density, since the expectation value of the others vanish. 
Given finite statistics, however, and an inexact estimate of the trace, it is possible that using the vacuum-subtracted 
value gives reduced statistical errors. Empirically, we find that for the strange axial form factor at Q 2 = 0, the 
two approaches give indistinguishable results. At larger momenta, however, uncertainties for the vacuum-subtracted 
quantities are noticeably smaller. 

To interpret our results, we require an understanding of the correlation functions given in Eqs. (I13l) - (|16[) in terms of 
the lowest one or two states that dominate at large times. This is accomplished by performing a spectral decomposition 
in the transfer matrix formalism, and for the nucleon two-point function given in Eq. (|12l) . we find 

G^(t, t ; = YJ 2 (l + |£) Z 2 n (q)e- E ^ . (19) 

Here Z n (p) is defined by (n,p, s\N a (0)\0) = Z n (p)u"(p), where N a (x) is a creation operator for the nucleon, \n,p, s) is 
its nth eigenstate (with momentum p and polarization s), and we have adopted a relativistic normalization convention 

for the states: (n',p',s , \n,p,s) — 2E n (p)L 3 d StS >d^,. The momentum dependence in Z n (p) arises because we will 
generally consider extended, rather than point-like, operators. We take n = 1 to label the ground-state proton, and 
for later convenience we collect together the coefficients, 

Cn{ q ) =2(l + ^Z 2 (q), (20) 

yielding 

G( 2 >(Mo;^ = £ C „(9>- £ '' (t - to) . (21) 

n 

Similarly, for a generic three-point function involving a current Jx(x, t), 

G^ 3) (M',*o;g) = rf YV^<A^(£,t)j x (a? ,i/)N a (o,t )) , (22) 

we find 

G^(t,t' > to;?1=E^(?) e_m " (t " t ' )e_Bm ® (t '^ o) ' ( 23 ) 

m,n 

where the coefficients j nm may generally be expressed in terms of suitable form factors, depending on the current 
Jx and the combination of gamma matrices Tx- The correlation functions in Eq. (|13p -(116 [) have been constructed 
such that in each case, the ground-state coefficient jn may be simply expressed in terms of a single form factor. In 
particular, for X = S, A, E, we have 



hi{q) = 1\l + j Z x {0)Z x {q)G x {Q') , (24) 

where G S X (Q 2 ) is the corresponding strange form factor of the nucleon. For the magnetic case, the appropriate 
expression is 

hi{q) = -^Z^ZMGIAQ 2 ) ■ (25) 
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FIG. 2. Fit of the nucleon correlation function for Q 2 — to a form that includes two forward-propagating states and one 
backward-propagating state (upper curve, red) and the same form with the coefficient of the backward-propagating set to zero 
but the other fit parameters held fixed (lower curve, blue). 



We work on a 24 3 x 64 anisotropic lattice, utilizing an ensemble of 863 gauge field configurations provided by the 
Hadron Spectrum Collaboration [50]. These were generated with two degenerate flavors in the sea. Our gauge and 
fermion actions are those defined in Appendix lAl with coupling (3 = 6/g 2 = 5.5 and bare anisotropy £o = 2.38. In [50l |. 
it was found that this value of £o together with v = 1 gives renormalized gauge and fermion anisotropies that are 
consistent with £ = 3. 

The spatial lattice spacing was determined from the Sommer scale [5l| with the parameter ro = 0.462(H)(4) fm, 
taken from [521 1 5 31] - yielding a s = 0.108(7) fm = 3at- For the light quarks, the mass parameter m® that appears in 
the action is m° = —0.4125. The corresponding pion mass is M„ — 416(36) MeV pSOj ]. 

Given the anisotropy, our lattice has a relatively short extent in time, which has influenced our choice of method. 
It has been conventional in lattice studies to extract the form factors by considering various ratios of the three- and 
two-point functions defined above. For example, at zero momentum transfer, one finds 



for large time separations. Instead, we have chosen to fit the three-point function that appears in the numerator of 
this ratio directly. The reasons are two-fold. First, this allows us to avoid contamination from backward-propagating 
states, which are problematic due to the short temporal extent of our lattice. At the same time, it allows us to 
explicitly take into account the contribution of (forward-propagating) excited states. 

To see why the direct approach avoids the problem of finite-time contamination, note that this contamination chiefly 
affects the nucleon correlator G^ (t,to;0) that appears in the denominator of Eq. (f26| . since it involves propagation 
for a time (t — to). As (t — to) exceeds L t /2, the correlator becomes progressively more contaminated by the negative- 
parity partner of the nucleon propagating backward through the lattice. In Figure [5J we show a plot of the nucleon 
correlator, together with a fit to a functional form that includes two forward-going states and one backward-going 
state. The lower curve shows the effect of dropping the term that corresponds to the latter; we use this for normalizing 
some of our results at zero momentum transfer when we plot them below. 

As a result of the contamination in the denominator, the ratio Rx begins to decrease precipitously at large times. 
It is important to note that although G x (t,t ,to;q) in the numerator also involves a nucleon propagating for time 
(t — to), the contamination there is a concern only insofar as it increases the statistical error by washing out the 
correlation we are attempting to measure. It remains unbiased since the current is inserted at t', while the negative- 
parity partner propagates across the opposite side of the lattice, from t to to, and can be expected to correlate little 
with the disconnected insertion. 



III. DETAILS OF THE CALCULATION 



Rx(t,t',t ;Q 2 = 0) 



G( 2 )(M ;0) 



^G x (Q 2 =0) 



(26) 
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FIG. 3. Nucleon energy-squared (in lattice units) as a function of momentum, together with a fit to the continuum dispersion 
relation. 



As described in the previous section, the correlation functions given in Eqs. p^)) - (IT5l) have been defined such that 
a single form factor enters the coefficient jn for each case, according to Eq. (l24l) . In terms of the coefficients c n 
extracted from the two-point function G^ 2 '(t, to', q), this becomes 



jn(q) = G s x (Q 2 )Jl (l + Cl (0) Cl (^ (27) 



. 2 I + E^q)) 

for X = S, E, A, where G s x (Q 2 ) is the corresponding strange form factor of the nucleon. The corresponding expression 
for G$(t,t',t ;q) is 



^-■^^f^mh^- (28) 

Our general strategy will be to fit the correlation functions G x (t, t' ,t ;q) to Eq. (1231) . taking into account both the 
ground state nucleon and a single excited state. We may then extract the nucleon form factors from jn with input 
from the two-point function. In principle, one could also obtain form factors of the first excited state from j'22, as well 
as transition form factors from ji2 and j2i- In practice, however, we expect these to absorb the contributions of still 
higher states and trust only the ground state form factors to be reliable. 



IV. RESULTS AND DISCUSSION 



A. Strange scalar form factor and fx s 



For all of the results presented in this section, the nucleon two-point function was fit in the range 10 < t/a t < 45, 
yielding a t M x — 0.198(2) for the ground-state nucleon mass. In Figure 02 we plot the nucleon energy (squared) as a 
function of momentum for the five smallest values of \p\ 2 available on our lattice, along with a fit to the continuum 
dispersion relation (atE) 2 = (a s |p|) 2 /£ 2 + (a t m) 2 . The parameter £, given by the inverse square root of the slope, 
provides a measure of the effective fermion anisotropy a s /at- We find £ = 3.25(11) from the fit, which may be 
compared to the values 2.979(28) and 3.045(35) obtained from the pion and rho dispersion relations, respectively, 
in [50|. The intercept (atm) 2 in Figure |3] is largely constrained by the point at \p\ 2 — and thus yields an identical 
value (and error) for the nucleon mass. 

In order to extract the form factors, the various three-point functions were fit to Eq. (|2"3"|) . taking into account the 
two lowest-lying states, with the separations (t — t')/dt and it 1 — to) /at varying independently in the range [10, 18]. It 



9 




FIG. 4. Subset of results for the scalar form factor at Q 2 — 0, where the current insertion is placed symmetrically between 
source and sink. The lower curve (red) shows a corresponding cross-section of the fit. The horizontal line (blue) indicates the 
resulting value of G%{Q 2 = 0)i a t = (N\ss\N)o for the ground-state nucleon. 
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FIG. 5. Strange scalar form factor as a function of momentum. 



follows that a total of 81 data points are included in the fit. A one- dimensional subset of these points for the scalar 
form factor at Q 2 = is shown in Figure [4] For the purpose of plotting, we have normalized our results by a fit 
to the two-point function. With this normalization, dominance of the ground state should manifest as a plateau at 
large times. We find G|(0)i a t = 2.84(49) for the form factor at zero momentum transfer, where the statistical error 
has been determined via a single-elimination jackknife applied to the full fitting procedure. In Figure [SJ we show the 
momentum dependence of the strange scalar form factor. 

We discuss the systematic uncertainties affecting these results below, including the delicate problem of relating 
the bare matrix element to the continuum. One practical consideration is the choice of fitting windows used in the 

(3) 

fits of the two- and three-point functions. In order to extract the form factors from G x (t, t , to; q), we must first 
determine the coefficients c n (q) and masses/energies E n (q) from a fit to G^(t,t ;q). Since we have access to a total 
of 863 x 64 — 55, 232 nucleon correlators, these tend to be very well-determined, as illustrated by Figure [2] The 
coefficients c„ are somewhat sensitive to the choice of fitting window, however, and since they multiply the form 
factor in Eq. (|27[) . this translates into a direct systematic error on the form factor, estimated to be about ten percent. 



10 



2 - 
1 





( At / a t)max= 16 

(At/a t ) mix =18 
(At/a,) max = 20 



10 12 14 16 



(At/a t ) n 



FIG. 6. Dependence of the extracted value of G|(Q 2 = 0)i a t = (iV|ss|iV)o on the range of time separations included in the fit. 




FIG. 7. The left diagram contributes to the renormalization of both the flavor singlet and flavor non-singlet mass operator, 
while the right appears only for the flavor singlet case. 



In contrast, we find that our results are relatively insensitive to the choice of window used in the fit of the three- 
point function. This is illustrated in Figure HI where we plot the extracted value of G S S (Q 2 = 0)i a t = (iV|ss|iV)o as a 
function of the smallest time separation included in the fit, for three different values of the maximum time separation. 
(See also the analogous plot for G S A (Q 2 = 0), Figure [T2l below.) We observe a stable plateau that extends to very 
early time separations but have nevertheless chosen a conservative lower bound, (i — t') > \0at and [t! — to) > 10a t , 
effectively eliminating systematics due to excited-state contamination of the three-point function, at the expense of 
increased statistical errors. 

Thus far, we have computed only the unrenormalized matrix element (JV|ss|iV)o on the lattice. Naively we 
can multiply by the subtracted bare strange quark mass fh s — mP s — m cr i t (determined in Appendix [Bj , to find 
(cs)iat = 504(91)(30) MeV, which in the continuum corresponds to the renormalization group invariant quantity 
a s = m s (N\ss\N). Here the second error reflects the uncertainty in the lattice scale, the first is statistical, and no 
other systematics have been taken into account. If we then divide by our measured value of the nucleon mass, a t Mj\r = 
0.198(2), the lattice spacing dependence drops out, yielding (fTs)ut — fh®(N\ss\N)o/MN — 0.46(9). This large value 
is in apparent disagreement with recent lattice determinations using staggered and chiral fermions [H, [13, IH, H(| ■ 
The source of this discrepancy, as first pointed out in [38] , is the explicit breaking of chiral symmetry in the Wilson 
action, which allows for mixing between singlet and non-singlet matrix elements even after tuning the quark masses 
— m® — ?7i cr it to zero. As a consequence, at finite lattice spacing the strange scalar matrix element receives 
contributions from both connected and disconnected diagrams involving the light quarks, such as those illustrated in 
Figure [7] 

As described in [54l.[55| in the context of quark mass renormalization, a natural approach for treating this problem is 
to separately consider the renormalization of flavor singlet and non-singlet contributions. In Appendix [Cl we employ 
this approach with the aid of the lattice Feynman-Hellmann theorem to rederive a result recently quoted in [37| | for 
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FIG. 8. Determination of a s — m 3 {N\ss\N) (in lattice units) from the discretized "kinetic term," as compared to same quantity 
evaluated by direct insertion of the scalar density multiplied by the subtracted bare quark mass. 



the renormalized matrix element, 

(N\ss\N) = 1 [{Z + 2Z s )(N\ss\N) Q + (Z - Z s )(N\uu + dd\N) ] + c(N\Tr[F 2 } \N) . (29) 

Here Zq and Z% are the flavor singlet and non-singlet renormalization constants for the scalar density, Tr[F ] is the 
gauge kinetic term, and c is a constant. The discussion in Appendix [Cl closely parallels the analysis of Bhattacharya, 
Gupta, Lee, Sharpe, and Wu [56|, who consider in detail operator mixing for Nf = 2 + 1 clover-improved Wilson 
fcrmions with = m u < m s , including all terms to O(a) and 0(am q ). This is a generalization of the classic on-shell 
0(a) improvement scheme of the ALPHA collaboration |57l. [58j. 

A self-consistent application of this approach demands an 0(a)-improved action with 2 + 1 dynamical flavors in the 
sea; such a calculation is under way (see conclusion) but beyond the scope of this paper. Nonetheless the discussion 
in Appendix \C\ is intended to clarify the source of the mixing problem. In accordance with [54l - [56| , it demonstrates 
that the singlet (Z™) and non-singlet (Z™) mass renormalization constants separately obey the reciprocal relations 
Z™ = 1/Zq and Z™ = \jZ% at zero quark mass and indicates why one expects the gluonic mixing (parameterized 
by c) to be small. In the approximation where the gluonic mixing is neglected (c — 0), correcting the dimensionless 
ratio /t s = m s (N\ss\N) /Mm only requires computation of the ratio Z s /Z . Consequently, we believe the value for 
the renormalization of the condensates can in principle be estimated following the prescription outlined in [Hil |55| ; 
by varying the valence and sea quark mass separately one can separate out the singlet and non-singlet contributions. 
Alternatively, one could determine this ratio by evaluating singlet and non-singlet matrix elements directly. Details 
of how best to compute the corrections are left for a future work. Suffice it to say that since the corrections due to 
mixing are large and negative (i.e., Z$/Zq > 1), we cannot rule out the possibility that the renormalized quantity 
(iV|as|JV) is consistent with zero within errors for the present calculation. 

We now consider an alternative method for determining <r s — m s (N\ss\N) by invoking the continuum equations of 
motion to replace m s ss by the quark "kinetic term," s^^D^s. On the lattice, the covariant derivative goes over to 
Dp = (Vyi + V*)/2, defined in terms of the covariant finite difference operators and V* given in Appendix [A"l 
Evaluating the matrix element of this operator gives us an alternative determination of a s , with lattice artifacts that 
are at least different and potentially less severe than those affecting the direct approach. In particular, this alternative 
obviates the need to separately consider quark mass and operator subtractions. Before turning to our results, we note 
that by splitting the lattice Wilson Dirac operator into three pieces, 

D = ^D f ,[U]+m l Q + W[U] (30) 

corresponding to the kinetic term, the bare mass term, and the Wilson term, respectively, and by invoking the exact 
lattice equation of motion, we may write this matrix element in two equivalent ways: 



- (N\s-ff,D^s\N) = (N\(m° s ss + sWs)\N) 



(31) 
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FIG. 9. Momentum dependence of a t m a G%(Q 2 ), as determined directly and from the kinetic term. 



From the expression on the right, we see that we have in effect subtracted the major shift due to TO cr it- 
alternative definition of the critical mass follows from imposing the condition, 



(N\(m° s ss + sWs)\N) = (m° s - rh clit )(N\ss\N) 



Indeed an 



(32) 



where m cr i t = —(N\sWs\N)/(N\ss\N). With this definition, at m° s — m cr it the flavor singlet term in Figure [7] is set 
to zero, suggesting that this scheme may suffer from smaller operator mixing than the direct approach. 

A practical question is how the statistical uncertainties in the two approaches compare. In Figure [8] we show our 
results for er s determined from the matrix element of the kinetic term. For comparison, we also include the data for 
(JV|ss|JV)o shown previously in Figure [4] but now rescaled by the bare subtracted quark mass (using the standard 
definition computed in Appendix IB"]) . If not for lattice artifacts, these two sets of results would correspond to the 
same continuum quantity. We find that the determination from the kinetic term does in fact suffer from much larger 
statistical errors, perhaps limiting the usefulness of the approach. A final judgement should await comparison of 
properly subtracted and renormalized results; such an investigation is underway For completeness, Figure [9] shows 
results for the scalar form factor as a function of momentum using the two approaches. The points with smaller error 
bars correspond to the data of Figure El rescaled by the bare subtracted quark mass. 

B. Strange axial form factor and As 

Results for the strange axial form factor are shown in Figure 1 101 here computed using the point-split current of 
Eq. (|18[). As was the case for G S S (Q 2 ), we note that our result (As)i a t = G A (0)i & t = —0.019(11) has not been 
renormalized and so may not be compared directly to experimental results!]] Despite the large errors, the data in 
Figure [10] seem to strongly favor a negative value for As, an observation that is in itself of phenomenological interest, 
given the present uncertainties in experimental determinations and the continued disagreement among some model 
calculations over the sign. In Figure fTTl we show the momentum dependence of G S A (Q 2 ). Because the determination of 
renormalization constants for our anisotropic lattice action is still pending, we present results for both the point-split 
and local axial currents. 

Finally, Figure fT2l shows how the extracted value of G S A (Q 2 = 0) would vary as a function of the range of time 
separations included in the fit. As described earlier, in computing all our results we have chosen a conservative range 
with the time separations (t — t')/at and {if — to) /at varying independently in the interval [10, 18]. The corresponding 
point in the figure is labeled by (At/a t ) m i n = 10 and (Af/a t ) max = 18. 



1 We also note that the preliminary results for the bare quantity G^(0)i at reported in |31| were computed using a point-split current 
involving gauge links rescaled by the bare anisotropy, £o = 2.38. Here we adopt a more conventional normalization for the current, for 
which Z\ ~ 1. 
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FIG. 10. Subset of results for the axial form factor at Q =0, where the current insertion is placed symmetrically between 
source and sink. The lower curve (red) shows a corresponding cross-section of the fit. The horizontal line (blue) indicates the 
resulting value of G^(Q 2 = 0)i a t = (As)i a t for the ground-state nucleon. 
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FIG. 11. Strange axial form factor as a function of momentum. 



C. Strange electric and magnetic form factors 

In Figure [13l we present our results for the strange quark contribution to the nucleon's electric and magnetic form 
factors, as a function of momentum. We have used the vector current defined in Eq. [TTl which is conserved for the 
Wilson action and therefore does not get renormalized. Note that since the strange quark does not contribute to the 
electric charge of the nucleon, G S E (Q 2 = 0) must vanish. This provides an additional check of our method, and we 
find G S E (Q 2 = 0) = —0.0016(20), consistent with zero as expected. More generally, all of our results for G S E (Q 2 ) and 
G S M (Q 2 ) appear to be roughly consistent with zero, implying that these quantities are rather small for Q 2 > 0.1 GeV 2 . 
Strictly speaking, we cannot set proper limits without extrapolating our results to the continuum and to the physical 
value of the light quark mass, but it is notable that the statistical errors are as much as an order of magnitude smaller 
than the corresponding experimental uncertainties (cf. [f| [H|). This suggests that measuring a nonzero value for 
G S E m(Q 2 ) m electron scattering experiments may be a challenging task indeed. 

In Table HI we summarize our results for the strange form factors of the nucleon, with the momentum transfer Q 2 
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FIG. 13. Strange electric and magnetic form factors as a function of momentum. 



given by 




where Mn is our lattice determination of the nucleon mass. The quoted errors for Q 2 reflect the uncertainties in 
Mat and the lattice scale. We again emphasize that the results in Table H] were determined with m u ,d unphysically 
heavy, corresponding to a pion mass of about 400 MeV, and that the tabulated values for G|(Q 2 ) and G S A (Q 2 ) arc 
unrenormalized. 



V. CONCLUSION 



In this work, we have described our first effort to compute disconnected contributions to nucleon form factors, 
focusing on the strange quark. Employing the Wilson gauge and fermion actions on an anisotropic lattice, we computed 
a large number of nucleon correlators and accurate unbiased estimates for the disconnected currents on each gauge 
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(L s /2^) 2 |gl 2 


Q 2 [GeV 2 ] 


Gs(Q 2 )iat 




Gi(Q 2 ) 1 ( :r al) 


Gi;(Q 2 )iat 










2.84(49) 


-0.019(11) 


-0.024(15) 


— 


— 


1 


0.22(3) 


1.59(34) 


-0.002(7) 


-0.008(10) 


0.000(3) 


-0.002(11) 


2 


0.43(5) 


1.20(28) 


-0.003(6) 


-0.006(9) 


-0.007(4) 


-0.005(8) 


3 


0.62(8) 


0.49(43) 


-0.001(11) 


-0.005(16) 


-0.003(7) 


-0.007(12) 


4 


0.81(10) 


-0.31(62) 


-0.012(15) 


-0.012(22) 


-0.009(11) 


-0.022(16) 



TABLE I. Summary of results for strange form factors of the nucleon. 



configuration. Wc nevertheless found results for the electromagnetic form factors that are consistent with zero and 
a result for As that is only marginally distinct from zero, suggesting that the physical values of these quantities are 
rather small. Such null results may be interpreted as limits — with the aforementioned caveats concerning systematics 
— and should also be useful for setting bounds on the disconnected contributions that are generally neglected in lattice 
determinations of nucleon form factors (or explicitly canceled by taking isovector combinations). To complete this 
program, it will of course be necessary to include disconnected contributions from light quarks as well. 

In the future, wc plan to build on the present investigation by introducing several improvements. First, we are 
making use of multiple ensembles of anisotropic lattices with 2+1 flavors in the sea [gcj, which will allow the strange 
quark to be treated fully self-consistently. These were generated with a Wilson fermion action that is stout-smeared [61j 
and 0(a)-improved [62j ]. both features that may be expected to improve the chiral properties of the action [63[ and 
thereby reduce the effect of flavor mixing discussed in Section flV Al Indeed, the fact that their action is clover- improved 
may explain why the authors of [35] found a value for (iV|ss|iV)o that is significantly smaller than ours (but still larger 
than determinations employing chiral or staggered fermions); one must also take multiplicative renormalization factors 
into account when comparing bare values obtained with different actions, but such factors are not expected to differ 
enough from unity to account for the discrepancy. These new ensembles also have a much longer extent in time (with 
volumes of 24 3 x 128 and larger), which will suppress contaminations from backward-propagating states and allow us 
to obtain a signal over a larger range of time separations, thus reducing statistical errors. 

Second, we are leveraging a powerful new adaptive multigrid (MG) algorithm for inverting the Wilson-clover Dirac 
operator that is allowing us to compute the disconnected diagrams for both strange and light quarks at very little 
additional cost [il, [65|. We are also taking advantage of clusters accelerated by graphics processing units (GPUs) 
using the QUDA library (6(| [63], which provides another substantial speedup. Work is underway to develop an MG 
implementation suitable for GPUs, in lieu of the Krylov solvers currently implemented in QUDA. We estimate that 
by combining these two improvements, we may be able to reduce the cost per Dirac inversion at light quark masses by 
up to two orders of magnitude as compared to standard solvers on traditional architectures. Finally, we are exploring 
additional methods for reducing the variance in estimates of the trace of disconnected currents, such as the multigrid 
subtraction method described in [45| . The net result of these improvements will be a significant reduction in both 
statistical and systematic errors. At the same time, the scheme outlined in Appendix [C] should allow us to correct 
for operator mixing in the determination of the strange scalar matrix element, yielding a reliable value and further 
elucidating the connection between results obtained with chiral and Wilson-like fermions. 
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Appendix A: The Wilson action for an anisotropic lattice 

In our calculation, we take the temporal lattice spacing a t to be finer than that in the three spatial directions, 
which share a common value a s . In the interacting theory, the anisotropy £ = a s /at renormalizes away from the bare 
value that appears in the action, which we denote by £o- Furthermore, the anisotropies appearing in the gauge and 
fermion actions may in principle renormalize differently. We follow [50] in denoting the bare gauge anisotropy by £o 
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while introducing a new parameter v such that the bare fermion anisotropy is given by £,a/v. We assume that the 
renormalization of the latter quantity is independent of quark mass, as found empirically in [Hoi. |69|. 
With these definitions, the Wilson gauge action on an anisotropic lattice is given by [70j 



Sn 



6 



EE 



x fj,— l \_fi<iy<4: 



(l-~B*U^(x)\ + $ (l-±BeU,*(x)\ 



(Al) 



in terms of the plaquette U^(x) = Tr [U ll (x)U^(x + a^fi)UUx + a v u)Ul{x)\ , where \i = 4 corresponds to the "time" 
direction. The Wilson fermion action, in turn, is given by [71j 



S w = «s E ^fc) 



/ 1 

a t m° q + -a 8 ]T ( - 7 *(V^ + V?) - -^V*V, 



<•, ( ^7 4 (v 4 +v:)-^v:v 4 



, (A2) 



where we have defined the covariant difference operators V M i/>(a;) = [U^(x)ip(x + a^fi) — tp(x)]/a^ and V* ll %l)(x) = 
[ip(x) — Ul(x — a fJi jj,)il>(x — a^jj,)]/ 'a^. Note that at occurs in Eq. (|A2[) only to cancel where it appears in the definition 
of V4, except in the dimensionless mass parameter (atmP q ). We can write the fermion action in a more familiar and 
explicit form by defining rescaled links, 



U il {x) = l €0 M 



f-U^x) for M = 1,2, 3 



and defining 



2 k 



[7 M (a;) for /.t = 4 . 



3// 



— = a t m" + — + 1 
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(A3) 



(A4) 



Thus we have Sw = a>l J2 X if>Dip(x), where 

1 1 4 

D ^ = _ 2 E [( : ~ Tm)^^)^^ + A) + (1 + 7/0^2(3 " A)^(a: - A) 



(A5) 



(For convenience, we have also redefined a^fi — > /t.) Note that we have kept the continuum normalization of the 
fermion field ip(x). 



Appendix B: Quark mass determination 

Due to the explicit breaking of chiral symmetry in the Wilson action, the naive quark mass m° that appears in 
Eq. (|A2I) is not protected from additive shifts under renormalization. In Section I IV A[ we require the subtracted bare 
mass, rhP s = mr s — m cl -it, where the critical mass m cr i t corresponds to the value of mo at which the physical quark 
mass vanishes. The naive mass for the strange quark is a parameter of the theory; it has been chosen such that the 
mass of the <j> meson calculated on the lattice reproduces the physical value [50| . To determine m cr i t , we utilize the 
dependence of the (partially quenched) pseudoscalar meson mass Mp on the valence quark mass: Mp = 2BmJ al to 
leading order in partially quenched chiral perturbation theory. The critical mass so defined depends implicitly on the 
fixed sea quark mass [72|, LZ3| ■ As advocated in [74[ , a reasonable approach is to determine m cr it for each available 
value of the light sea quark mass 77i| ea and then extrapolate to the physical point where mp a = m Ut d- (Alternatively, 
one could fit all available data for Mp to a functional form that incorporates the dependence on both the sea and 
valence quark masses [55[.) Consistent with the other results presented in this work, because only a single value of 
?7ij oa is available, we do not perform this final extrapolation in the light sea quark mass. 

Figure [TJ] illustrates our determination of the critical mass. For each of seven valence quark masses, we evaluated 
pseudoscalar correlators on an ensemble of 216 gauge configurations. The corresponding meson masses were deter- 
mined from a single-cosh fit in the range 20 < tj at < 44. Upon performing a linear extrapolation in the quark mass, 
we find atm CY it = —0.42116(24), where the statistical error has been estimated via jackknife. Given the naive mass 
that was input for the strange quark, a t m° = —0.38922, we find atm° s = 0.03194(24) for the subtracted bare strange 
quark mass. 
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Appendix C: Flavor mixing for Wilson quarks 

In principle to extract continuum quantities, one must take the lattice spacing (i.e., the bare coupling) to zero 
holding renormalized parameters fixed, even for renormalization group invariant quantities such as ratios of masses. 
However, better estimates can often be found at finite lattice spacing by "renormalizing" the bare lattice quantities. A 
particularly interesting and difficult quantity for Wilson fermions is the continuum parameter for the Higgs coupling 
to the strange quark content of the nucleon, 

m s (N\ss\N) d , 

St. = \j n = m s — log[Mjv] . (CI) 

The expression on the right is an identity based on the Feynman-Hellmann theorem with the partial derivative taken 
with respect to the renormalized strange quark mass, holding the renormalized light quark mass and scale fixed. For 
Wilson quarks on the lattice, mixing with the light quark condensates in the nucleon can produce a large contribution 
to Jts, as pointed out by Michael, McNeile, and Hepburn (38|. To understand this, let us consider a mass-independent 
renormalization scheme on the lattice. We note that very similar methods are employed in Section IV. B of [56], despite 
some differences in the choice of the lattice renormalization scheme. For convenience, in this Appendix, we adopt a 
convention where lattice mass parameters are dimensionless. Also we will consider only the more physically relevant 
case of 2 + 1 flavors, even though our present results are calculated on 2 flavor gauge configurations. There is however 
still non-zero mixing in the 2 flavor case as illustrated in Figure [7] to lowest order. 

For Wilson quarks there are two important issues not present for a chiral formulation. First, we have additive mass 
renormalization, which requires that a constant be subtracted from the bare quark mass. Second, the disconnected 
diagram for the mass insertion operator tptp = ipLipR + 4 , r4 , l does not vanish even when the subtracted quark masses 
vanish. We begin by rewriting the mass term in the lattice Lagrangian in terms of singlet, m° s — (2mf + mj;)/3, and 
non-singlet, m° NS — (rnf — m®)/y/3, masses: 

C m = m®(uu + dd) + m s ss = rrigipip + rn^vs^MsV' , (C2) 

where for simplicity we take degenerate light quarks = rrP u = mj. The singlet and non-singlet masses renormalizc 
differently because of the lack of chiral symmetry. The renormalization scheme we choose is 

m s = Z™(g )(m° s - m crit {go))/a , 
m NS = Z , 8 n (g )m° NS /a, 

A = A°(g )/a. (C3) 

The first two equations define the renormalized masses, while the last defines some renormalized scale which, for 
simplicity, we will take to depend on the lattice spacing and some function of only the bare coupling. There are many 
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possible choices for A, such as the rho mass or pion decay constant, though it will not enter into our final results so 
we leave it unspecified. We have also introduced the lattice spacing, a, to convert quantities to physical units. This 
needs to be set by comparing some lattice measurement with a physical value. One choice would be to set the lattice 
spacing using the Sommer scale, (5l| . as 



-rt yS /r l Q at i9o) 



(C4) 



where rg at is the dimensionless value measured on each lattice ensemble and rQ hys is some reference value in physical 
units. As with A, the exact choice of definition for a is irrelevant for the present discussion. Here we have made 
these quantities independent of the masses. One could systematically improve on this by adding extra terms with 
increasing powers of the subtracted masses, but, for simplicity, we will not include these. 

The mass renormalizations above can be expressed in terms of the light and strange quarks themselves, 



mi =m s + m NS 



/VZ = \ [(2Z™ + Z™)m° + (Z™ - Z?)m°] /a , 



m s = m s - 2m NS /Vz = ^ [(Zj? + 2Z™)m° s + 2(Z™ 



»?] A*. 



(C5) 



where = m° — m cr i t . We are now ready to find the expression for renormalized condensates. Due to the vector 



Ward identity the non-singlet operator renormalizes as 



(C6) 



where Zg = l/Z™, but the singlet piece is unconstrained by vector current conservation. However we can uniquely 
determine the renormalization of the condensates by evaluating the Feynman-Hcllmann theorem, 
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(C7) 



(where Mjy = aMjsr is the nucleon mass in lattice units) in terms of bare lattice parameters. The partial derivative is 
expanded as 
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leading to the expression 
dm® 



h 



dm. 



dm s 



(^|ss|A^)o 



dmP s d dm® d dg 2 d 



dm s dm , dm s dm ' dm s dg 



dm. 



■(N\uu + dd\N) Q 



dm s 



(N\g 2 Q S g \N) a 



a ; M 



(C8) 



(C9) 



where S g is the gauge action and (.)o is an unrenormalized lattice matrix element. 

Finally, using the renormalization scheme in (IC3[) . we evaluate the coefficients in this expression by use of the 
implicit function theorem, inverting the Jacobian matrix 
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From Eq. (|C5I) the determinant is then given by 
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Matrix elements of the inverse of the Jacobian are thus given by 
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Thus we identify Zq = 1/Z™ and Z$ — l/Z™ to obtain a form similar to Eq. (|29l) in the main text, 

frs = 3^ [ (Z + 2Z s ){N\ss\N)q + (Zq - Z 8 )(N\uu + dd\N) ] . (C13) 

It is interesting that here the relations Zi = l/Z™ did not involve the use of Ward identities, contrary to standard 
derivations. Also note that there is no (N\ g^S g \N)Q contribution to this order in the renormalization scheme. However 
as emphasized in (56| . it may be important to include additional 0(am s ) corrections which will cause dg^ 2 / dm s to 
no longer vanish and will induce operator mixing with (N\Sq\N)q. Since on dimensional and RG grounds this term is 
0(am s g^), it should be relatively small. The mixing with the light valence quarks is substantial and with the estimate 
for Z$/Zo > 1, it will tend to cancel the contribution from (iV|ss|iV)o found for the bare amplitude. To see this in 
more detail, note that renormalizing Jt s also requires finding the renormalized strange quark mass, 

™ s = \ [(Z n + 2Z™)m° s + 2(Z™ - Z™)m°] /a . (C14) 

Consequently, the dimensionless ratio fa. only depends on the ratio Z™/Z™ = Z s /Z , which can be computed by a 
procedure described in detail in [H H in 

the context of quark masses. It is convenient to rewrite the expression as 
f Ts = ^ + 2A(m°- m°)/[3(l + A)] _ A{jv|( ^ + ^ _ 2Sa)]N)o/3 ] , (C 15) 

M N 

in terms of an operator mixing parameter A = Z$/ 'Zq — 1. The mixing term is a pure non-singlet operator, A-ipXgip/ V3- 
The disconnected contribution to ^Ag^/vS vanishes like 0(ami) in the chiral limit, but the valence contribution 
remains large at current lattice spacings, resulting in a correction of the same order of magnitude as the bare matrix 
element given in the text: a- 1 fh s (N\ss\N) ~ 504(91)(30) MeV. 
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